GENERIC CONTROLLABILITY OF 3D SWIMMERS IN A PERFECT FLUID 
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Abstract. We address the problem of controlling a dynamical system governing the motion of a 3D weighted 
shape changing body swimming in a perfect fluid. The rigid displacement of the swimmer results from the exchange 
of momentum between prescribed shape changes and the flow, the total impulse of the fluid-swimmer system being 
constant for all times. We prove the following tracking results: (i) Synchronized swimming: Maybe up to an 
arbitrarily small change of its density, any swimmer can approximately follow any given trajectory while, in addition, 
undergoing approximately any given shape changes. In this statement, the control consists in arbitrarily small 
superimposed deformations; (ii) Freestyle swimming: Maybe up to an arbitrarily small change of its density, any 
swimmer can approximately tracks any given trajectory by combining suitably at most five basic movements that 
can be generically chosen (no macro shape changes are prescribed in this statement). 
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1. Introduction. 

1.1. Context. Researches on bio- inspired locomotion in fluid have now a long history. Focus- 
ing on the area of Mathematical Physics, the modeling leads to a system of PDEs (governing the 
fluid flow) coupled with a system of ODEs (driving the rigid motion of the immersed body). The 
first difficulty mathematicians came up against was to prove the well-posedness of such systems. 
This task was carried out in [16] (where the fluid is assumed to be viscous and governed by Navier 
Stokes equations), in [T3] (for an inviscid fluid with potential flow) and in [S] (for low Reynolds 
numbers swimmers, the flow being governed by the stationary Stokes equations). 

After the well-posedness of the fluid-swimmer dynamics were established, the following step 
was to investigate its controllability. On this topic, still few theoretical results are available: In [2], 
the authors prove that a 3D three-sphere mechanism, swimming along a straight line in a viscous 
fluid is controllable. In [4|, we prove that a generic 2D example of shape changing body swimming 
in a potential flow can track approximately any given trajectory. 

Some authors are rather interested in describing the dynamics of swimming in terms of Geo- 
metric Mechanics (within the general framework presented for instance in |12| ) . We refer to [7] and 
the very recent paper [8] for references in this area. 

In this article, we consider a 3D shape changing body swimming in a potential flow. Under 
some symmetry assumptions (the swimmer is alone in the fluid and the fluid-swimmer system fills 
the whole space) we prove generic controllability results, generalizing and improving what has been 
obtained for a particular 2D model in [?]. 

1.2. Modeling. 

Kinematics. We assume that the swimmer is the only immersed body in the fluid and that 
the fluid-swimmer system fills the whole space, identified with R"^. Two frames are required in the 
modeling, the first one 2; := (Ei,E2,E3) is fixed and Galilean and the second one e := (ei,e2,e3) 
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is moving with its origin lying at any time at the center of mass of the swimming body. At 
any moment, there exist a rotation matrix R G S0(3) and a vector r g R'^ such that, if X := 
{Xi,X2tX^Y and x := (a;i,X2,X3)* are the coordinates of a same vector in respectively £ and e, 
then the equality X = Rx + r holds. The matrix R is supposed to give also the orientation of the 
swimmer. The rigid displacement of the swimmer, on a time interval [0, T] (T > 0), is thoroughly 
described by the functions t : [Q,T] ^ R{t) e S0(3) and t : [Q,T] ^ r{t) g R^, which are the 
unknowns of our problem. Denoting their time derivatives by R and r, we can define the linear 
velocity v := (wi, W2, 1^3)* G R-'^ and angular velocity vector ft := (f^i, il2, f^a)* G R"^ (both in e) 
by respectively v := i?*r and fl := R^R, where for every vector u := (ui,W2,W3)* G R'^, u is 
the unique skew-symmetric matrix satisfying ua; := u x a; for every x G R'^. Vectors of R^ will 
sometimes be decomposed in the form f := (f^,f^)* G R'^ x R'^. For every f (f^,f^) G R^ and 
g := (g\g^) € R^ we can define f Tk-g := (f^ x P,f^ x - gi x f^)* G R^. 

Shape Changes. Unless otherwise indicated, from now on all of the quantities will be ex- 
pressed in the moving frame c. In our modeling, the domains occupied by the swimmer are images 
of the closed unit ball B by diffeomorphisms, isotopic the the identity, and tending to the 
identity at infinity, i.e. having the form Id -I- ■)? where i? belongs to Z)g(R^) (definitions of the 
function spaces are given in the appendix. Section With these settings, the shape changes 
over a time interval [0, T] can be simply prescribed by means of an absolutely continuous function 
t G [0,T] !->■ G £'q(R'^). Then, denoting 0t = Id -I- i^f, the domain occupied by the swimmer at 
the time t > is the closed, bounded, connected set Bt := 0t{B) (do not forget that we are working 
in the frame e). We still require some notation: the unit ball's boundary is E := dB, Et := 6*4(1]) 
stands for the body-fluid interface, nt is the unit normal vector to Et directed toward the interior 
of Bt and the fluid fills the exterior open set J-j := R"^ \ Bt- 

So-called self-propelled constraints are necessary to ensure that the deformations result from 
the work of internal forces (they avoid for instance translations to be considered as allowable shape 
changes). Let a function g G C^{B)^ be given. The density of the deformed body at the instant t, 
denoted by gt G C°{Bt), is defined by: 

gt{x):=g{e^\x))/Jt{0^\x)), {x e Bt, t > 0), (1.1) 

where Jt := |detV0t| = det V0t (we can drop the absolute values here because 0t{x) x as 
I|2^||r3 +00 and deiVOt{x) ^ for all x G R"^ and t > 0). The self-propelled constraints read: 

/ gt{x)xdx = and / gtix)dtetiO^^{x)) x xdx ^ (t > 0). (1.2a) 

The former identity means that, as already mentioned before, the center of mass of the swimmer 
lies at any time at the origin of the moving frame. The latter relation tells us that the angular 
momentum (in e) has to remain constant as the swimmer undergoes shape changes. Equivalent 
formulations can be obtained up to a change of variables: 

/ g{x)et{x)dx = and / g{x)dtet{x) x 6>t(x)dx = (t > 0). (1.2b) 

J B J B 

The Flow. The fluid is assumed to be inviscid and incompressible. We denote hj gf > 
its constant density. The flow is governed by Eulcr equations. According to Helmholtz's third 
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Fig. 1.1. Kinematics of the model: The Galilean frame £ := (Ej)i<j<3 and the moving frame c := (ej)i<j<3 
with Bj = REj (R G SO{3)). Quantities are usually expressed in the moving frame. The domain of the body is Bt 
at the time t and Bt is the image of the unit ball B by a diffeomorphism Ot- The open set Tt '■= R.'^ \ Bt is the 
domain of the fluid. The center of mass of the body is denoted r (in £j, v := /?*r is its translational velocity (in () 
and O its angular velocity. 

theorem, if the flow is irrotational at the initial time, it remains irrotational for aU times. In this 
case, the Eulerian velocity is equal at any time to the gradient of a potential function. According to 
Kirchhoff's law, the potential can be decomposed into a linear combination of elementary potentials, 
each one connected to a degree of freedom of the system (they consist here in the 6 degrees of freedom 
of the rigid motion of the body plus those connecting to the deformations). These ideas have been 
thoroughly described in a series of papers [3 [131 [Ml |4] , to which we refer for further explanations. 

The elementary potentials "0? j (* = Ij • ■ • 1 6) corresponding to the rigid motion of the swimmer 
are harmonic in J-t, tend to zero as infinity and satisfy the Neumann boundary conditions dn'fpl = 
(e^ X x)-r\.t {i = 1, 2, 3) and dnipl = e.i-^-nt {i = 4, 5, 6) on Sf. They are well defined in the weighted 
Sobolev space W^{J^t) (defined in the appendix. Section [Al see also [3j for details). The overall 
potential connecting to the rigid displacement is Vt := J2i=i ^ii^t + Si=4^i-3V't (* > 0). On the 
other hand, the elementary potential ipt associated to the shape changes, harmonic as well in J^t, 
satisfies the boundary condition dn(pt = Wt ■ rit on Sj {i = 1, . . . ,n), where Wt(a;) := dt0t{O^^{x)) 
{x e R^). Like the functions tpl {i = 1, 2, 3), (ft belongs to W'^{Tt) for all t > 0. 

Dynamics. The modeling of moving rigid (or shape changing) bodies in an ideal fluid classi- 
cally involves the notion of mass matrices. The mass of the body is m :~ Jg gdx and its inertia 
tensor at the time f > is defined by: 



We introduce MJ,'(<) diag(I(f), mid) (a 6 x 6 symmetric bloc diagonal matrix), Mj(f) (a 6 x 6 
symmetric matrix as well) whose entries read: 



I{t) := / gtlWxW^sld - X (g) x] dx 




Q[\\Ot{xWR,ld - Ot{x) ® Ot{x)] dx. 



(1.3) 



JBt 




(1.4) 
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and we denote WF'{t) := MJJ(t) + M^(t). We also need the 6x1 column vector N(t), homogeneous 
to a momentum, whose elements read: 

Qf I Vipl-Viftdx, {l<i<6). (1.5) 

If we neglect the buoyancy force, it has been proved in a series of papers (we refer for instance to 
the already mentioned articles [71 [T3] or 0]) that the swimming motion is governed by the equation: 



(^^) = -M'-(t)-iN(i), it > 0). (1.6a) 

At this point, we can identified gt (or more simply g since they are linked through relation (|l.ip ) 
as a parameter and the control as being the function t G [0,T] i-> G Dq(R'^). Notice that the 
dcpendance of the dynamics in the control is strongly nonlinear. Indeed 'dt describes the shape of 
the body and hence also the domain of the fluid in which are set the PDEs of the potential functions 
involved in the expressions of the mass matrices M''(t) and N(t). 

To determine the rigid motion. Equation (|1.6a[) has to be supplemented with the ODE: 



d_ /RX _ fRQ, 
dt\r J [rv 



{t > 0), (1.6b) 



together with Cauchy data for i?(0) and r(0). Remark that the dynamics does not depend on g 
and gf independently but only on the relative density g/gf. So we can assume, without loss of 
generality, that gf ^ 1 in the sequel. 



1.3. Main results. The first result ensures the well posedness of System ()1.6p and the conti- 
nuity of the input-output mapping: 

Proposition 1.1. For any T > 0, any g g C^{B)^ , any absolutely continuous function 
t G [0,T] t— > G £'o(R^) (respectively of class , p = l,...,-t-oo or analytic) and any inital 
data (i?(0),r(0)) G S0(3) x R^, System (fOl) admits a unique solution t G [0,T] ^ {R{t),r{t)) G 
S0(3) X R'^ (in the sense of Caratheodory) absolutely continuous on [0, T] (respectively of class 
or analytic). 

Lett G [0,r] H- G Dl{I{?) (for) ^l,...,+oo) be a sequence of controls in AC {[0,T], Dl(R^)) 
(see Section\^for a definition of this space) which converges in this space to a function t G [0, T] i— > 
■dt e I?o(R-^)- Let a pair (i?o,ro) G S0(3) x R^ be given and denote t G [0,r] ^ _{R{t),r{t)) G 
S0(3) X R3 the solution in AC([0, T], S0(3) x R^) to System (HH]) with control and Cauchy 
data (i?o,ro). Then, the unique solution {R^,r^) to System (|1.6p with control "d^ and Cauchy data 
(i?o,ro) converges in AC([0, T], S0(3) x R^) to {R,r) as j +oo. 

We denote by M(3) the Banach space of the 3x3 matrices endowed wit any matrix norm. The 
main result of this article addresses the controllability of System (jl.6p : 

Theorem 1.2. (Synchronized Swimming) Assume that the following data are given: (i) A 
function g in C'^{B)^ (the target density of the swimmer); (ii) A C^ function t G [0,T] i— ?> i^t G 
Z)q(R'^) (the target shape changes) such that the pair {g,'&) satisfies the self-propelled constraints 
([TS)) .- (Hi) A C^ functiontG [0,T] ^ {R{t),f(t)) G SO(3)xR3 (the target trajectory to be followed). 
Then, for any £ > 0, there exists a function g G C^{B)^ (the actual density of the swimmer) and a 
function t G [0,T] i— >^ i9t G D^iTi?) (the actual shape changes that can be chosen of class C^ for any 
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p = 1, . . . , +00 or even analytic) such that the pair {g, i?) satisfies (|1.2[) . Wq — £'||c"(S) < £7 '^o = ^^O; 

dj, = iJj, and suptg[o.T] (ll^t " llci(R.-'')^ + ~ ^(i)llM(3) + II^W - i"(i)llRO < £ where the 

junction t G [0,r] t-^ {R(t),Y[t)) G S0(3) x R'^ is the unique solution to system (jl.6p wit/i initial 
data (i?(0),r(0)) = (i?(0),f(0)) and control 

This theorem tells us that, maybe up to an arbitrarily small change of its density, any weighted 
3D body undergoing approximately any prescribed shape changes can approximately track by swim- 
ming any given trajectory. It may seem surprising that the shape changes, which are supposed to 
be the control of our problem, can also be prescribed. Actually, be aware that they are only ap- 
proximately prescribed. We are going to show precisely that arbitrarily small superimposed shape 
changes suffice for controlling the swimming motion. This result improves what has been done in 
the article [J] for a particular 2D model. 

When no macro shape changes are prescribed we have: 

Theorem 1.3. (Freestyle Swimming) Assume that the following data are given: (i) A pair 
(^, 1?) S C'^{B)~^ X £)q(R^) satisfying (|1.2|) (the target density and shape at rest) (ii) A function 
t G [0,r] 1-^ {R{t),v{t)) G S0(3) X (the target trajectory). Then, for any e > there exists a 
pair (p, I?) e C°(i?)+ X Dq(R^) (the actual density and shape at rest) such that (i) Jg gOdx ~ 
(where := Id + d) (ii) \\q — q\\c°{b) + ll^^ ~ ^1IdJ(r3) < £ '^'^^ (Hi) for almost any 5-uplet 
(Vi, . . . , V5) e (Co^(R3)3)5 satisfying gV.dx = 0, gO xV.dx = and gV, x V^- dx = 
(i,j = l,...,5j, there exists a function t G [0, T] t— > s{t) G R^ (that can be chosen of class 
for any p = 1, . . . , +00 or even analytic) such that, using "dt := + X)i=i Si{tyVi G D}^{Y{?) as 
control in the dynamics (jl.6p . we get sup^gjQ (||_R(t) — i?(i)||]vi(3) + ||r(t) — r(t)||f{^3) < e where the 
function t G [0,T] n> {R{t),r{t)) G S0(3) x R^ is the unique solution to ODEs (|1.6p with initial 
data (i?(0),r(0)) ^ (i?(0), f (0)). 

Differently stated, wc claim in this Theorem that any weighted 3D body (maybe up to an 
arbitrarily small change of its density) is able to swim by means of allowable deformations (i.e. 
satisfying the self-propelled constraints) obtained as a suitable combination of pretty much any 
given five basic movements. 

The proofs rely on the following main ideas: First, we shall identify a set of parameters necessary 
to thoroughly characterize a swimmer and its way of swimming (these parameters are its density, 
its shape and a finite number of basic movements, satisfying the self-propelled constraints (jl.2[) ). 
Any set of such parameters will be termed a swimmer configuration (denoted SC in short). Then, 
the set of all of the SC will be shown to be an (infinite dimensional) analytic connected embedded 
submanifold of a Banach space. 

The second step of the reasoning will consist in proving that the swimmer's ability to track any 
given trajectory (while undergoing any given shape changes) is related to the vanishing of some 
analytic functions depending on the SC. These functions arc connected to the determinant of some 
vector fields and their Lie brackets (we will use here some classical result of Geometric Control 
Theory). Eventually, by direct calculation, we will prove that at least one swimmer (corresponding 
to one particular SC) has this ability. An elementary property of analytic functions will eventually 
allow us to conclude that almost any SC (or equivalently any swimmer) has this property. 

Remark 1.4. The authors conjecture that in both Theorem \1.2\ and Theorem ] 1.31 the actual 
density g can be chosen equal to the target density g. At this point however and although it is very 
unlikely, it can not be excluded that all of the swimmers with a particular density might be unable 
to swim. This issue also appeared in 
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1.4. Outline of the paper. The next Section is dedicated to the notion of swimmer con- 
figuration (definition and properties) . In Section [3] we show that the mass matrices are analytic 
functions in the SC (seen as a variable) and in Section|4]we will restate the control problem in order 
to fit with the general framework of Geometric Control Theory. A particular case of swimmer will 
be shown to be controllable. In Section[5]the proof of the main results will be performed. Section|6] 
contains some words of conclusion. Many technical results are gathered in the appendix to avoid 
overloading the rest of the paper. 

2. Swimmer Configuration. A swimmer configuration is a set of parameters characterizing 
swimmers whose deformations consist in a combination of a finite number of basic movements. 

Definition 2.1. For any positive integer n, we denote C(n) the subset ofC^{B)~^ x _Dg(R-^) x 
(Cq(R"^)'^)" consisting of all of the triplets c := {g,'d,V) such that, denoting :~ Id + ■& and 
V :~ (Vi, . . . , V„), the following conditions hold (i) the set {V^j^ -e^, I <i < n, k ~ 1,2,3} is a 
free family in Cq{B) (ii) every pair (V, V) of elements o/{6?, Vi, . . . ,V„} satisfies gV dx = 
and Jg gV X Y' dx = 0. 

We call swimmer configuration (SC in short) any element c of C{n). 

By definition, £)g(R^) is open in C^i^^Y (see appendix. Section [^. We deduce that for any 
c G C{n), the set {s := (si, . . . , s„)* £ R" : d + Ylll=i ^i^i ^ ^o(^^)} open as well in R" and 
we denote 5(c) its connected component containing s = 0. 

Definition 2.2. For any positive integer n, we call extended swimmer configuration (ESC in 
short) any pair c := (c, s) such that c £ C{n) and s G S{c). We denote Cxin) the set of all of these 
pairs. 

Restatement of the problem in terms of SC and ESC. Pick a SC in C{n) (for some 
integer n). The characteristics of the corresponding swimmer can be deduced from c as follows: 
If c is equal to (g, i9, V) with V := (Vi, . . . , V„), the shape of the swimmer at rest is B := 0(B) 
where := Id + "d. When swimming, it can occupy the domains Be '■= 0s{B) for all s G S{c) (c := 
(c, s) G C{n) is hence an ESC), where 0s ■= ^d + d + SiYi. Still for any s G S{c), its density is 
the function g^ G C°(Sc)+ defined by g^ix) ~ g{0-'^{x))/Jsi0j'^{x)) with dot VQ^. Notice 

that within this construction, the shape changes on a time interval [0, T] (T > 0) are merely given 
through an absolutely continuous function t : [0,T] i-^ s{t) G 5(c). If i G [0,T] i-^ s{t) G R" 
stands for its time derivative, the Lagrangian velocity at a point a; of i? is X)i=i Si(t)Yi{x) while 
the Eulerian velocity at a point x £ Be is where w*(a;) := Yi{0J^ (x)). Due to 

assumption (ii) of Definition 12.11 the self-propelled constraints (|1.2p are automatically satisfied. 

The harmonic elementary potential functions of the fluid corresponding to the rigid motions 
depend only on the ESC. Therefore, they will be denoted in the sequel f/'c instead of ipl- The 
same remark holds for the inertia tensor I{t) and the mass matrices M''(t), MJJ(i) and My:(t) whose 
notation is turned into 1(c), M''(c), MJJ(c) and M[j(c) respectively. The elementary potential 
connected to the shape changes can be decomposed into Y^'^=i ^iVc- this sum, each potential 
function ipl. is harmonic in := R'^ \ Be and satisfies on Sc := dBc the Neumann boundary 
conditions dnfl. = • Uc , Uc being the unit normal to Ec directed toward the interior of Be- 
Introducing the mass matrix N(c), whose elements are gf Vtpc ■ S/(pldx (1 < « < 6, 1 < j < n) 
(recall that gf can be chosen equal to 1), the dynamics (|1.6ap can now be rewritten in the form: 

= -M'Xc)~iN(c)s, (t > 0). (2.1) 
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Let us focus on the properties of C{n) and Cxin). 

Theorem 2.3. For any positive integer n, the set C{n) is an analytic connected embedded 
submanifold of C°{B) x CqHR^)^ x {C^(R^f)" of codimension N := 3(n + 2)(n + l)/2. 

The definition and the main properties of analytic functions in Banach spaces are summarized 
in the article [TTj . 

Proof For any c := (£>,i?, V) G C"{B) x Co1(R3)3 x {C^{R^)^y\ denote Vq := Id + i9 and 
V := (Vi,...,V„). Then, define for k ^ 0,1,..., n, the functions : C"(B) x C^iK^f x 

l3^3^n 



(Ci(R3)3)» ^ R3(n+i-fc) by Afc(c) := ( J^gVkdx, J^gVk x V^+i da;, . . . , x V„ da 

Every function is analytic and we draw the same conclusion for A := (Aq, . . . , A„)* : C°(i?) x 
Co^R^)^ X (CqHR^)^)" R^ (N := 3(n + 2)(n + l)/2). In order to prove that dcA{c) (the 
differential of A at the point c) is onto for any c G C{n), assume that there exist (n + 2)(n + l)/2 
vectors aj G R'^ (0 < i < j < n) such that: 



1=0 



a, ■ {d,K{c),c^) = 0, Vc'* G C°{B) x Co'(R^)^ x (Coi(R^))^ 



where ai := {ol\,ol\ 



(2.2) 



Ci(R^)^ X (Ci(R3)3)" with V(J 
(12.21). wc obtain that: 



G R3(«+i-») {j = 0, ...,n) and c'' := {q'\{)^,V'') G C°{B) x 



Id + and V'' := (V^*, . . . , Vjj). Reorganizing the terms in 



?''[x:«^Vfe+ ^ «i.(v, xv,) 

0<i<j<n 



fe=0 



di- 



rt k~l n 

fc=0 j=Q j=k+l 



dx = 0. 



Since this identity has to be satisfied for any (p'', i^'', V'') G C"{B) x C^CR^f x (CoHR^))^ we 
deduce that, for every k = 0, . . . ,n: 



E 

p=0 



aP ■ V 



p + E 

o<j<i<Ti 



a^ -(V, xVj)=0 and 



fe-i 

E 

J=0 



-E ' 

J=fc+1 



X V, 



(2.3) 



Multiplying the latter equality by g and integrating over B, we get that = (fc = 0, . . . ,n). 
The first equality in p.3p is now just a linear combination of the second ones (for k = 0, . . . ,n), 
so we can drop it. Taking into account Hypothesis (ii) of Definition 12.11 latter identity in p.3p 
with k = leads to ckq = for every j = 1, . . . ,n. There are no more terms involving Vo in the 
other equations and invoking again Hypothesis (ii) we eventually get aj=0 for 1 < i < j < n. 
So, equality p.2p entails that cti = for alH = 0, . . . , n and the mapping dcA{c) is onto for all 
c G C{n). 

The linear space X = Ker9cA(c) is closed since A is analytic. Let Y be an algebraic supplement 
of X in C"{B) X C^in^f x (Co^(R3)3)", and denote by Hy the hncar projection onto Y along X. A 
crucial observation is that the linear space Y is isomorphic to R^ and hence it is finite dimensional 
and closed in C°{B) x Co^(R^)^ x {C^(R^)^Y\ Define the analytic mapping / : AT x F -> R^ by 
f{x,y) = A(c + X + y). The mapping dyf{0,0) = 9cA(c) o Hy being onto, the implicit function 
theorem (analytic version in Banach spaces, see jl7j ) asserts that there exist an open neighborhood 
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Oi of ill X, an open neighborhood O2 of in Y, and an analytic mapping g : Oi Y such that 
g{0) = and, for every {x, y) in O1XO2, the two following assertions are equivalent: (i) f{x, y) = 
(or, in other words, c + x + y belongs to C(n)), and (ii) y = g{x). The analytic mapping g provides 
a local parameterization of C{n) in a neighborhood of c. 

In order to prove that C{n) is path-connected, consider two elements := {Q\d\V'^) and 
ct := (e*,t9t,V*) of C(n) and denote := Id + := (Vj,...,Vt) and 6>* := Id + z?*, 

V* := (Vf, . . . , V* ). According to Definition US £I^(R3) is open and connected. It entails that 
it is always possible to find a continuous, piecewise linear path t : [0, 1] i-> G Dq{Y{?) such that 
"01=0 = and •dt=i = i^''"- We introduce = < <i < . . . < t/c = 1, a subdivision of the interval 
[0,1] such that t ^ dt \s linear on every subinterval {j = 0, — 1) and we denote 

Ot ■= Id + := ^t=tj , 0-' Id + ^^-^ (j = 0, . . . , fc). Let us introduce as well the continuous 

functions t £ [0, 1] 1-^ gt := tgi + {l-t)Q^ £ C°{B) and t G [0, 1] H- u(t) := - £»t(9tda;/ /g ptdx £ 
R'^. The set Utg[o.T]0t(-B) + u(t) being compact, it is contained in a large ball il. We introduce 
r2' an even larger ball containing ft and a cut-off function x defined in R'^ such that < x < 1; 
X = 1 in and x = in R^ \ fj'. The derivative ii of u exists everywhere on ]0, T[ excepted maybe 
at the points ti, . . . The flow associated with the Caratheodory's solutions of the Cauchy 

problem: Xt{x) = u{t)x{Xt{x)), {t > 0), Xt=o{x) = a; is well defined (see [TTJ Theorem lA, page 
57]). Moreover, for every fixed t G [0, 1], the mapping x G R'^ 1— > Xt{x) G R^ is C°°. Consider now 
the mappings t G [0, 1] ^^ i?t_:= Xt o Ot - Id_and := Id -I- -dt. If x G R^ \ Cl', 0t{x) = Otix) 
for ah t G [0,r] and if x e B then -dtix) = i?t(a;) + u{t) and ©((cc) = Otix) + u(t). Notice that 
-dt e £'o(R^) for ah t G [0,T] and /g gtix)et{x)dx = for ah t G [0,T]. Since Co^R^)^ is an infinite 
dimensional Banach space, it is always possible to find by induction Wi, W2, . . . , W„ in Co(R^)'^ 
such that (i) both families {Wi|g • e^, . . . , W„|^ • efe, V||^ • efc, . . . , |g • efc, k 1,2,3} and 
{Wils-Gfe, . . . , W„|s-efe, Vfl^-efc, . . . , V* l^-efc, k = 1,2,3} are free in C|J(R^) and (n) for any pair 
of elements V, V'( both picked in the same family, g^Ydx = 0, Jg g^Vdx = 0, Jg g^O^ x Vdx = 
0, /g gi& X Vdx = (for ah j = 1, . . . , fc), g'^V x V'dx = and /g gW x V'dx = 0. Define 

the function t G [0, 1] Vj G C^CR^f by Vj := (1 - 2t)V| + 2tW, ff < t < 1/2 and Vj := 
(2 - 2t)W, + (2t - l)Vi ff 1/2 < t < 1 and denote Vt := (V^, . . . , V^") G {C^(R^)'^)". Eventually, 
a continuous function linking to is given by t G [0,1] 1— S- Cf G C{n) with Cf := (f?^, i?^, V3j/2) if 
<f < 1/3, Q := (^)3^-l,^?3t-l,Vl/2) if 1/3 < t < 2/3 and q :== (^)^ , i?*, V3t/2-i/2) if 2/3 < t < 1. □ 

We omit the proof of the following corollary, similar to that of the theorem above: 

Corollary 2.4. For any positive integer n, the setCx{n) is an analytic connected embedded 
suhmanifold of C"{B) x C^{Ii^)^ x (Co^R^)^)" x R" of codimension N := 3(n -|- 2)(n + l)/2. 

We denote tt the projection of C{n) onto C^{B) x Dq{T{?) defined by 7r(c) — {g,'d) for all 
c := {g,-d,V) G C{n). The proof of the following corollary is a straightforward consequence of 
arguments already used in the proof of Theorem 12.31 

Corollary 2.5. For any positive integer n and for any (gj-d) G 7r(C(n)), the sectiomT~^{{{g,'d)}) 
is an embedded connected analytic submanifold of {g} x x (Cq^RP)^)" (identified with (CqCRP)^)" ) 
of codimension 3n{n + 3)/2. 

3. Sensitivity Analysis of the Mass Matrices. For any positive integers k and p, we 
denote M{k,p) the vector space of the matrices of size k x p (or simply M(fc) when fc = p). 

Theorem 3.1. For any positive integer n, the mappings c G Cxin) 1— > M[^(c) G M(6) and 
c G Cx{n) M> N(c) G M(6,n) are analytic. 

The method followed in this proof is inspired from [6] . The result already appeared in [13j , in 
a slightly different form though. Due to its crucial importance for our purpose, we recall here the 
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main ideas. 

Let us begin with a preliminary lemma of which the statement requires introducing some 
material. Thus, we denote F := \ B (recall that B is the unit ball, S := dB and n is the 
unit normal to S directed toward the interior of B). For all ^ S DqCEL^), we set S := Id + ^, 
;= ^{B), J^^ := S{F), := S'(E) and stands for the unit normal to directed toward 
the interior of B^. We denote q := (C, W), with W := (W\ W^) c {C^{R^ f)'^, the elements of 
Q := DlCR^) X (Co1(R3)3)2 and := [i = 1,2). Finally, for every q e Q, we define: 

<P{q) := f ^,pl^ix)-Vi^l{x)dx, (3.1) 

where, for every i = 1,2, the function ipl^ g W^{T^) (recall that the function spaces are defined 
in Section is solution to the Laplace equation —Aipl^ = in J-'j with Neumann boundary data 
dnipq " • on E^. The solution has to be understood in the weak sense, namely: 

/ S/4>\{x)-Vip{x)dx^ I {wi-n^){x)ip{x)da, £ VF^J^f). (3.2) 

Lemma 3.2. The mapping q G Q ^ ^(q) G R is analytic. 

Proof. We pull back relation (|3.2p onto the domain F using the difFcomorphism S. We get: 
/ V*';(x-)A^(x)-V(^oS)(x)dx= / {W\x)-n^{S{x)MS{x))J^{x)da, \/^eW\T^), 

J F JT. 

where if* := Vq o ^, dct(VH), := {\/S*VE:)-'^J^ and := |[(VS-i)*n||R3 (usually 

referred to as the tangential Jacobian). In (|3.2p . if we specialize the test function to have the 
form ip (j) o with G M^^(F), we obtain Jp^^{x)A^{x) ■ V(/)(x) dx = /j, &q(x)0(a;) dcr for 
all (/) G W^{F), where fe^ := (W* • n^(S'))Jf (i = 1,2). We now claim that the mapping ^ G 
D^CR^) i-^ - Id G E^{n, M(3)) is analytic. Indeed, the mappings ^ G DlCR^) i-> VS^VS - Id G 
i;°(f],M(3)), A G E^{n,M{3)) ^ (Id+A)-i-Id G E^{n,M{3)) and^ G £'ii(R^) ^ J^-l G Co"(R3) 
are analytic. Reasoning the same way, we can show that the mapping ^ G Dq(R^) m> G C*'(E) is 
analytic as well (notice that for all ^ G Z)q(R'^), the function (VS^^yn never vanishes on E). It is 
more complicated to prove that ^ G -Dq(R'^) i-> linos' G C°(E)^ is analytic, so we refer to [T^ for the 
details. This last result entails the analyticity of q G Q &q G C°(E). Then, denoting by W^{Fy 
the dual space to W^{F), we consider the mapping T : (q, m) G Q x W^{F) {A^,u, ■) - (6q, ■) G 
PFi(F)', where (A^, := VM(x)A^(a;)-V(/)(a;) dx and := /j. &jj(a;)<?!)(a;) da {(j)eW\F)). 

We wish now to apply the implicit function theorem (analytic version in Banach spaces, as stated 
in |17| ) to the analytic function F. Observe, though, that we are only interested in the regularity 
result and not in the existence and uniqueness. Indeed, we already know that for every q G Q, 
there exists a unique function G W^{F) such that r{q,'F^) = 0. The function if'q is equal to 
tpl^o S where lAq is the unique solution to the well posed variational problem p.2p . The partial 
derivative 9„P(q, tf'q) is defined by: 

(a„r(q,if^),w,0) = / Vii(a;)A^(x) • V</)(a;)dx, V0 G I^^(i^). (3.3) 
Jf 

For all ^ G -Dq(R'^), the matrix A^ is uniformly elliptic in R'^ (there exists > such that 
X*A^{x)X > a^||X||^3 for all X G R^ and all x G R^). Wc deduce that the right hand side of 
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can be chosen as the scalar product in W^{F) and hence that 9tiP(q, tf'q) is a continuous 
isomorphism from W^{F) onto its dual space according to the Riesz representation theorem. The 
implicit function theorem applies and asserts that the mapping q G Q i— > •f'q £ W^{F) is analytic. 

To conclude the proof, it remains only to observe that the function <?(q) introduced in p.ip can 
be rewritten, upon a change of variables as ^(q) = JpW'F^{x)A^{x) ■ V'F^lx) dx, which is indeed 
analytic as a composition of analytic functions. □ 
We can now give the proof of Theorem 13.11 

Proof. Recall that the elements of the matrix Mj(c) are defined in (|1.4p and those of N(c) in 
(II. 5p . For any c (c, s) G Cx{n), where c := {g, z?, V), we apply the lemma with ^ :— ^+Y^^=i ^i'^i 
and W\W2 e {e, x S, e^, i = 1,2,3} to get that the mapping c G Cx{n) ^ M'^(c) G M(6) 
is analytic. To prove the analyticity of the elements of N(c), wc apply the lemma again with 
^ ■= + YT^=lSi^^^ G {e, X E, e„ i = 1,2,3} and G {Vi,...,V„}. Eventually, the 
analyticity of the elements of MJJ(c) is straightforward after rewriting the inertia tensor 1(c) defined 
in (|1.3p in the form (upon a change of variables) 1(c) = £)[||S'||^3ld ~ S ^ S] dx, still with 
S:=ld + C and ^ 7? + ^Li ^^V,. □ 

4. Control. 

4.1. Controllable SC. Let us fix c G C{n) (for some positive integer n) and recall that 5(c) 
is the connected open subspacc of R" such that (c, s) G Cx{n)- Introducing (fi, . . . , f„) an ordered 
orthonormal basis of R", we can seek the function t G [0,T] i-^ s{t) G 5(c) as the solution of the 
ODE s{t) = X;r=i Ai(i)fi where the functions A, : i G [0,T] t-^ X,{t) G R are the new controls, and 
rewrite once more the dynamics p.l|) as: 




-M'^(c,s)-iN(c, s)f, 
f,; 



{t > 0). (4.1) 



It is worth remarking that in this form, s is no more the control but a variable which is meant to be 
controlled and c G C{n) is a parameter of the dynamics. Considering (|4.ip . we are quite naturally 
led to introduce, for all c G Cxin), the vector fields X,(c) := -M''(c)-iN(c)f, G R*', Y,;(c) := 
(Xi(c),X2(c),f,)* G ridS0(3) X R3 X R" (we have used here the notation X, := iXj,Xfy G 
R3 xR3) and Z*(i?,s) := 7eflY,(c) G rfl,S0(3) x R^ x R" wh ere TZ r := diag(i?, i?. Id) G S0(6 + n) 
is a bloc diagonal matrix. The dynamics (|4.ip and the ODE (jl.6b[) can be gathered into a unique 
differential system: 



d 
di 



= ^A.(OZ^(i?,s), it>0). (4.2) 



For every i = 1, . . . , n, the function (i?, r, s) G S0(3) x R^ x 5(c) ^ Zi{R, s) G TflS0(3) x R'"^ x R" 
can be seen as an analytic vector field (constant in r) on the analytic connected manifold A4{c) := 
S0(3) x R'^ X 5(c). We denote any element (i?, r, s) G A4{c) and we define Z{c) as the family of 
vector fields (Z*)i<i<„ on M{c). 

Lemma 4.1. Let c be a SC fixed in C{n) (n a positive integer). If there exists C G A4(c) such 
that dimLie^.E(c) = 6 + n, then the orbit of Z{c) through any C, G M.{c) is equal to the whole 
manifold M.{c). 

Proof. Rashevsky Chow Theorem (see [T]) applies: If Lie^Z(c) = TqM{c) for all C G M{c) 
(or more precisely, for all (i?, s) G S0(3) x 5(c) since does not depend on r) then the orbit of 
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Z{c) through any point of M[c) is equal to the whole manifold. Let us compute the Lie bracket 
[Z*(i?,s), ZJ(i?,s)] for l<i,j <n and {R,s) e S0(3) x <S(c). We get: 



[Zi{R,s),ZiiR,s)]=.nR 



\ 




.Xi)(c) 

(XI X x2 - X] X xf )(c) I + 7^fl I (a,,x^ - a.^.xf )(c) I . (4.3) 



By induction, we can similarly prove that the Lie brackets of any order at any point C S Ai{c) 
have the same general form, namely the matrix TZr multiplied by an element of T'(id,o,s)-^(c). We 
deduce that the dimension of the Lie algebra at any point of A4{c) depends only on s. According 
to the Orbit Theorem (see [T]), the dimension of the Lie algebra is constant along any orbit. But 
according to the particular form of the vector fields Z* (whose last n components form a basis of 
R"), the projection of any orbit on S{c) turns out to be the whole set <S(c) (or, in other words, 
for any s € S{c) and for any orbit, there is a point of the orbit for which the last component is s). 
Assume now that dimLie^.Z(c) = 6 + n at some particular point := {R* ,r* , s*) G Ai{c). Then, 
according to the Orbit Theorem, for any s £ 5(c), there exists at least one point (i?s, r^, s) G A4{c) 
such that dimLie(fl^^r„,s)'2^(c) = 6 + n. But since the dimension of the Lie algebra does not depend 
on the variables R and r, we conclude that dimLie^Z(c) = 6 + n for all ( € Ai{c). □ 

Definition 4.2. We say that c, a SC in C'{n) (for some integer n) is controllable if there 
exists C G 7W(c) such that dimLiefZ(c) = 6 + n. 

It is obvious that for a SC to be controllable, the integer n has to be larger or equal to 2. The 
following result is quite classical (a proof can be found in [3]): 

Proposition 4.3. Let c G C{n) (for some integer n) be controllable (with the usual notation 
c := (q, 1?, V), V := (Vi, . . . , V„) and := "d + Ylll=i ^i^i for every s G S{c) ). Then for any given 
continuous function t G [0,r] t-^ {R(t),r{t), s(t)) G S0(3) x x S{c) and for any £ > 0, there 
exist n functions Ai : [0, T] R (i ~ 1, . . . ,n) such that: 

1. SUPtg[o j-j {\\m - R{t)\\M{3) + l|r(i) — r(t)||R.'i + \\i}s{t) — '^s(t) II ci(R3)3) < 

2. i?(T) = R{T), r{T) = f(T) and s(T) = s{T); 

where t G [0,T] (R(t),r(t), s{t)) G M{c) is the unique solution to the ODE (|4.2|) with Cauchy 
data R{0) = R{0) G S0(3), r(0) = f (0) G R^ s(0) = s(0) G 5(c). 

Let us mention some other quite elementary properties that will be used later on: 

Proposition 4.4. 

1. Ifc:^ ((?,i9, V) G C{n) (n > 2) is a controllable SC with V := (Vi, . . . , V„) G {C^(Il^)^)" 
then any c+ := ig,^,V+) G C{n+1) such that V+ := (Vi, . . . , V„, V„+i) G (CoHR^)^)"^^ 
(for some V„_|-i G C(P (Rp)^ ) is a controllable SC as well. 

2. If c := (p,i9, V) G C{n) (n > 2) is a controllable SC, then for any i?* G + s fVi, s G 
5(c)} the element c* :~ (p, V) belongs to C{n) and is a controllable SC as well. 

3. If c := (p, i9, V) G C{n) (n > 2) is a controllable SC, then all of the controllable SC 
c* := (p, 1?, V*) in C{n) (V* G {Cq(R.^)^)" ) form an open dense subset of the section 
"""^({(Pi i^)}) (for the induced topology). 

4. If there exists a SC in C{n) for some n > 2 then, for any k > n, all of the controllable SC 
in C{k) form an open dense subset of C{k) (for the induced topology). 

Proof. The two first assertions are obvious so let us address directly the third point. Denote 
Sk {k positive integer) the set of all of the vectors fields on A4{c) obtained as Lie brackets of order 
lower or equal to k from elements of Z{c). Then, consider the determinants of all of the different 
families of 6 + n elements of S^. as analytic functions in the variable V (the other variables g, i? and 
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s = being fixed). Since c is controllable, there exist at least one k and one family of 6 + n elements 
in £k whose determinant is nonzero. According to Corollary 12.51 and basic properties of analytic 
functions (see |17|). the determinant may vanish only in a closed subset with empty interior of the 
section TT~^{{{g,-d)}) (for the induced topology). The proof of the last point is similar. □ 

4.2. Building a controllable SC. In this subsection, we are interested in computing the Lie 
brackets of first order [Z* (i?, s), Zl{R, s)] at (i?, s) = (Id, 0), for a particular SC. 

General computations. Starting from identity (|4.3p and focusing on the second term in the 
right hand side, we have, for all c G C{n) (n > 2) , all s G S{c) and all i, j = 1, . . . , n: 



a,.X,(c) - 9,^X,(c) = M'^(c)-^ (9,^.M''(c)X,(c) - 9,,M''(c)X,(c)) 



+ {ds^c^ - dsMc)fj) ■ (4.4) 



From the decomposition M''(c) := M^c) + M^.(c), we deduce that: 



9,^M'^(c)X,(c) - dsM^{c)X,{c) = 

idsMl{c)X,{c) - dsMl{c)X,{c)) + {dsM}{c)X,{c) - dsM}{c)X,{c)). (4.5) 

We can easily compute the first term in the right hand side when s = 0. Thus, for all i, j = 1, . . . , n, 
we have: 

[5,^M^(c)X,(c) - a,,M^(c)X,(c)],=o = V, 

where i]i := 2 g{x) Vi(x) ■ 0{x) dx. 

Rigid shell's deformation. We consider now shape changes that reduce to rigid displace- 
ments on the swimmer's boundary S (but obviously not inside the body for the self-propelled 
constraints not to be violated). The idea of using such deformations stemmed from the reading 
of the article [S]. However, rigid deformations of the shell can not be considered within the mod- 
eling described in Section [5] (because rigid deformations do not fit with the general form of the 
diffeomorphisms described in Section [2]). Let us explain how to deal with this difficulty. 

Let c := {g,'&,V) £ C(6) be given such that V := (Vi, . . . , Vg) with V,|s(a;) := e, x {x + d{x)) 
{i = 1,2,3) and Vi|s(x) := ei_3 (i = 4,5,6) (Proposition IC.ll in the appendix guarantees the 
existence of such vector fields satisfying furthermore rji = for alH = 1, . . . , 6). Thus, the shape 
changes we are considering read 0s := Id -I- -I- X]i=i ^iVi (s G S{c)). Let us define also '■= 
i?s(Id + + Ts where Rs := exp(siei) exp(s2e2) exp(s3e3) G S0(3) and Ts := Z^Li ^i+^^i G R-^- 
Thus, 0s is a diffeomorphism which can be achieved within our modeling while Ss is a true rigid 
deformation. 

In order to determine the expressions of the terms 9s.My-(c)|s=o (* = 1, ■ • • , 6) arising in the 
computation of the Lie brackets, we have to compare (using the notation of Lemma l3.2|) 9s.«?(q^)|s=o 
and dsi'P{c^s)\s=o (* = I7 • • • i 6) where: 

• q.' (el, Wi) with ■■= -Id, Wi := (W^^i, W^'^) and Wl'\ W^^^ ^ ^ ^ ^ 
1, 2, 3} (settings corresponding to a true rigid deformation of the shell); 

• q' := i^l^^s) withe? := 0.-Id = ^ + ^^V,, := {Wf\Wf^) SindWf\Wl^ G 
{ei X 0s, ei, i = 1,2,3} (settings allowed in our modeling). 
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Id 




Xi(c)|s=o 



Id 




X,(c)|,=o, (4.6) 



Remark that ds^Osls^a ~ dsi'^s\s=o on S for all i = 1, ... ,6, so the deformations are tangent at 
s = 0. It entails that 9siqsls=o = dsiq^\s=o- Since, in addition, q^^o = ^8=0 ^^'^ 9si^{q'^)\s=a = 
(5qf?(q^),5s,qs)|s=o for fc = 1,2 and i = 1,...,6, we deduce that ds,${ql)\s=o = 9s,<?(q^)|s=o for 
alH = 1, . . . , 6. 

We apply the diffeomorphism Ss := Rs{ld + 'd) + Ts to B to obtain the domain of the deformed 
swimmer. We denote := {ld + -d){B) (which can be seen as the shape of the swimmer at rest, i.e. 
when s = 0), := 96*, J"* \ and Bs := E:,{B), \ S„ E, dBs. We seek the 

potential ip^, defined and harmonic in Ts in the form V'c(^) = V'K^sl^ ^ '''s)) where the function 
ipl defined on has to be determined. It is obvious that ipl is harmonic in J^* and we have only 
to determine the boundary conditions on E*. For all j/ G E*, we denote x := Rgy + Tg G E^. The 
relation n I s,(i?s2/+Ts) = Rsn\Y,o{y) entails that Vi/i^(y)-n| so (y) = Vipl{Rsy+Ts)-n\^^{Rsy+rs) = 
S/tpKx) ■n{x), and this quantity has to be equal to (ei x x)-n{x). We deduce after some elementary 
algebra that dnfpliy) = {Rlei x y) ■ n(j/) + (i?*ei x RIts) • n(y). Proceeding similarly and with 
obvious notation, we obtain more generally that: 



dn^liy) 
and 9„Vi*(y) 



(i?*e, x y) ■ n{y) 
Rlet-3 ■ n{y), 



(i?*e, X i?*r,) • n(y). 



(i = 1,2,3), 
(i = 4,5,6). 



Denoting :— Mj(c)|s=o, we get the relations: 
and then, differentiating with respect to s^: 

a,,M;^(c,s)|,=o 



"Id Ts' 




'Rs ■ 




Rl 


■ 




'Id 0" 


Id 




R, 







Rl 




-T, Id 



"e, 0' 






e. 0" 


e, 






e, 





e,_3' 






0" 


— 









ei-3 



(i = 1,2,3), 
(z = 4,5,6). 



(4.7a) 
(4.7b) 



The reasoning is quite similar for the entries of the matrix N(c). We only have to be care- 
ful that the vector fields W^^^ can actually not depend on s (once more, this is due to our 
modeling). The elements of the matrix 9s N(c)|s=ofi ~ 9s;N(c)|s^ofj read (still with the nota- 
tion of Lemma 1321) ds^<P{q^,)\s=o - ds,<Piqi)\s=o = 1,...,6) with this time ql := (Cs,W]), 
Cs := - Id, Wl := (Wi, W'-2), G {e^ x ~„ e^, k = 1,2,3} and := x (Id + i?) 

(i = 1,2,3) or = e,_3 {i = 4,5,6). The only difference with the elements of M^(c) 

being that W'^ docs not depend on s for i = 1,2,3, we wish to reuse the preceding calcula- 
tions. Invoking the chain rule, we have to subtract to dsj'P{ql)\s=o — 9si^(qi)|s=o the quantity 
{^^^2<l>{ql) , dsjWl'^)\s=o — (9w2^(qi), 9^. W^'^)|s=o which is quite easy to determine because ^ is 
linear in the variable W^. Thus, this last expression is merely equal to «?(qj;.j) (1 < i,j < 6) with 
q,,, := (^?,>V,,j), W^,J := (W^Wf A, G {e,. x (Id + i9), e^, k = 1,2,3} and 



W"^ - •= 



['(e,; X Gj) X (Id + 1?) if 1 
(e,; X e,) 



<i,j < 3, 
if 4 < i < 6 , 1 < j < 3, 
if 1 < j < 3,4 < j < 6, 



if 4 < i, j < 6. 
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We eventually obtain that: 



a,^N(c)|,=of. ~ 5,.N(c)|,=of, = d,M}{c)\s=oi^ - a,.M^(c)|,=of,- - N*(f, *f,), (4-8) 

where := N(c)|5=o = and * f, := (f^ x fj,fi x f| - fj x iff with the notation f, = 
(fi,f2)* eR3 xR3. 

Specifying the density and shape. The expression of the 6x6 symmetric added mass 
matrix depends only on the domain J^* or equivalently on B*. As stated in Proposition ID . 1 1 in 
the appendix, this matrix is positive definite if we choose i? in such a way that has no axis of 
symmetry. It entails that, up to a change of frame e at the initial time, Mj can be assumed to be 
diagonal with positive eigenvalues (j = 1, . . . , 6). On the other hand, denoting by 8(3)+ the set 
of the 3x3 symmetric matrices that are positive definite, we can quite easily prove that for any 
(which means for any ■& G Dl{R^)), the mapping g g C°(S*)+ ^ /^^ £i(||a;||R3ld-a:(g)a;)da; G S(3)+ 
is onto. We deduce that for any {Ii,l2,h) G such that Ii > for i ~ 1,2,3, there exists 
g G C°(i3)+ such that the inertia tensor I(c)|s=o is diagonal, equal to diag(/i, /2, /s). In this 
case, the matrix M^(c)|s=o is diagonal as well, equal to diag(/i, /2, /s, to, m, to, ). We deduce that 
the vector fields Xi(c)|s=o read —fii/{Ii + iJ.i)fi if i — 1,2,3 and ~^i/{m + ^i)fi if i = 4,5,6. 
Summarizing prt]) . (|i75|) (rccaU that 77^ = for ah i = 1,...,6), and P^ . we 

obtain the following expressions for the Lie brackets at (i?, s) = (ld,0): 



'63' 



,yl 731 ^ hhit^i - M2 + M3) + MlM3(A - ^2 + ^3) 

^ ^' (M1+/i)(A^2+/2)(M3+/3) 

hhifJ-i - Ai2 - M3) + M2M3(A - -?^2 - -^3) 



[Z^, z^] 



[zi Z^] 



\K. Z^] 



(Ml +/i)(M2 +/2)(M3 +^3) 

/2to(m4 - /ie) 



(M2 


+ /2)(Ai4 + 


m){iJ.e - 


\-m) 






- M4) 




(M3 


+ /3)(Ai4 + 


to)(jU5 - 


\-m) 






- M4) 




(M3 


+ h)in4 + 


w)(m5 - 


\-m) 



Since, on the other hand, when (i?, s) ~ (Id, 0) we also have 
-Miej/(Mj + Ii 

I if 1 = 1,2,3, and Zl = \ -^iie,/{fi, + m) \ ifi = 4,5,6. 
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we deduce that a sufficient condition ensuring that dimLie(id,o.o)'2^(c) = 12 is that 



[hhi-fJ-i - M2 + A's) + yui/i2(-/i - h + h)] [^^3(^1 - M2 + A's) + yUiA'a^i ~ I2 + h)] 

[hhifH - - Ma) + ^^2^^?,{h - h- h)] [hm{^n - ^iQ)llm^{^iz - ^4)^] 7^ 0. (4.9) 

According to [TOl pages 152-155], if we specialize now to be an clhpsoid, the length of the axes 
can be chosen in such a way that (i) it has no axis of symmetry (and hence /i^ > for i = 1, . . . , 6), 
(ii) /i4 ^ /is and (iii) /i4 ^ /ig. Since Ii > (i = 1,2,3) and obviously m > 0, the condition (|4.9p 
reduces to: 

[/i/2(-/ii - M2 + A'a) + MiM2(--fi - ^2 + 4)] [^1^3(^1 - M2 + A*3) + MiA*3(-^i - -^2 + 4)] 

[^2/3(^1 - M2 - M3) + M2M3(/i -h- h)] 7^ 0. (4.10) 

As already mentioned, it is always possible to achieve any triplet of positive numbers {Ii, 12,13) 
with a suitable choice of density, so whatever the values of Hi [i = 1, 2, 3) are, it is always possible to 
find g € C°{B)+ such that holds. It entails, according to the forth point of Proposition 14.41 

Proposition 4.5. For any integer n > 5, the set of all the controllable SC is an open dense 
subset in C{n). 

Notice that n = 5 in this Proposition (instead of n — 6), because we did not use the vector 
field in the computation of the Lie brackets. 

5. Proofs of the Main Results. 

Proof of Proposition II. li Let a control function d be given in AC{[0,T], DqCR?)) and 
denote := Id + i?. With the notation of Lemma 13. 2[ at any time t the entries of the matrix 
M^(i) read ^'(q) with q := (i?,>V) and W {W^,W^), W* e {e^ x 6>t, e^, i = 1,2,3}. We 
deduce that t e [0,T] ^ M^(i) G M(3) is in ^C([0, T], M(3)). To get the expression of the 
elements of the vector N(i) we only have to modify which has to be equal to dt^f It entails 
that t e [0,T] 1^- N(t) e R*" is in L^([0,T])'^. It is quite easy to verify that the inertia tensor 
I{t) is in AC{[0,T],M{3)). Existence and uniqueness of solutions is now straightforward because 
t G [0, T] ^ M'-(f)-iN(i) G R^ is in AC{[0, T], R^). 

Let us address the stability result. With the same notation as in the statement of Proposi- 
tion [TTll denote (f2-',v-')* the left hand side of identity (|1.6a|) with control i?-' and (J7,v)* when 
the control is ^. As j +00, it is clear that (fl-'^v^)* (f2,v)* in ^^([O^r], R^). Then, in- 
tegrating (jl.6bp between and t for any < i < T, we get the estimate \\R{t) — i?-'(0llM(3) ^ 
/p ll-R(s) — i?-'(s)||M(3)||i^(s)|jR3 + ||f^"'(s) — f2(s)||R3ds. Applying Gronwall inequality, we conclude 

that — > R in C([0, T], M(3)) as j +00 and we use again the ODE to prove that W R in 
L^([0,T]). Next, it is then easy to obtain the convergence of r^ to r and to conclude the proof. 

Proof of Theorems 11.21 and 11.31 We shall focus on the proof of Theorem 11.21 because it 
will contain the proof of Theorem 11.31 For any integer n, we shall use the notation |jc||c(„) := 
McoiB) + ll^llc„i(R3)3 + ELi l|V.||c,i(R3)3 for aU c G C°{B) x C^iR^)' x {C^iR^fT with, as 
usual, c := {g, V) and V := (Vi, . . . , V„). 

Let e > and the functions g G C°{B), t G [0, T] ^ e DKR^) and t G [0, T] ^ {R(t),r{t)) G 
SO (3) x R'^ be given as in the statement of the theorem. 

Set now g^ := g, := dt=o and V} := dt^t=o & C'g(R'^)'^. According to the self-propelled 
constraints (|1.2p . it is always possible to find four elements Vj (j = 2,..., 5) in Cg(R'^)'^ such 

15 



that the element := V^) be a SC which belongs to C(5), with := (V}, . . . , V^). Then, 

Proposition l4.5l guarantees that for any (5 > it is possible to find a SC controllable in C(5), denoted 
by := (g^-dW^) where := {Yl...,Vl), such that \\c^ - ci||c{5) < S/2. Since t ^ dt^t is 
continuous on the compact set [0,r], it is uniformly continuous. For any > 0, there exists > 
such that Wdtdt — dt'0t'\\c^{R^)^ < ^ providing that \t — t'\ < 5^. We then divide the time interval 
[0, T] into = ti < t2 < ■ . ■ < tk = T such that \tj+i ~ tj\ < 6^ for j = 1, . . . , - 1. For any 
ti < t < i2, we have the estimate: 

\\dt - + {t- tl)V})||cl(R3)3 < W^t - i^' + it~ tl)V})|bl(R3)3 

+ 11^' - ^ic„l(R3)3 +{t- h)\\Vl - Vl||cl(R3)3. 

On the one hande, we have, for all t € [^1,^2], \\^t — {^^ + {t ^ ^i)VJ)||ci(r3)3 < i^\t — ti\. On the 
other hand, still for ti < i < ^2 and if we assume that 6^ < 1, we get: — '&^\\cl,{'R^)'^ + ~ 
^i)ll^i ~ "^1 llc^(R'')3 ^ <5/2- We introduce now g'^ :— g^ and i?^ := ^t=t2- It is always possible to 
supplement := dt^t2 with vector fields V| (j = 2, . . . , 5) in such a way that := {g'^^-d^, V^) 
be in C(5) with the obvious notation := (V^, . . . , Vg). 

We define g^ := and := + (^2 — For any ti < t < t2, Proposition 14.41 guarantees 

that the SC := {g^,^^ + ^ ^1)^17^^) is controllable. In particular, for t = t2, there exists 
an integer k and a family of 11 vector fields in Sk (the set of all the Lie brackets of order lower 
or equal to k) such that the determinant of the family is nonzero. But this determinant can be 
thought of as an analytic function in V^. The set TT~^{{{g^ j-d^)}) being an analytic connected 
submanifold of (Co(R'^)^)^ (see Corollarv l2.5p . the determinant is nonzero everywhere on this set 
but maybe in a closed subset of empty interior (for the induced topology) . Therefore, it is possible 
to find G (Ci(R3)3)5 such that \\c^-c'^\\e{5) < (3/2 + iy{t2 ~ h)) + 5/4, and := {g'^,^^,V'^) is 
controllable. 

By induction, we can build and cJ {j ^1,2,..., k) such that ||c^' - || < 6/2 + Y/i=2 -^/S* + 
v{ti — ti_i) < (5 + vT. We choose 5 and v in such a way that 5 + vT < e/2 and we define 
t : [0,r] I— ??t 6 _Dq(R'^) and t G [0,T] n> ct G C(5) as continuous, piecewise affinc functions by 
dt := d^^(t-tj)N\ and q := (^, i9t, with ~g = g\ V* := ift G [tj, t^+i] (j = 1, ■ . ■ , k-1). 

Notice that for any t G [0,r], (i) W-fft — '!?t||c^(R3)3 < e/2 and (ii) q is controllable. 

Definition 14.21 and Proposition 14.31 ensure that, on every interval [tj,tj^i] (j = 1, . . . , fc — 1), 
there exist five functions : [tj,tj+i] i-)- R (i = 1,...,5) such that the solution {Rj,rj,s^) : 
[tj,tj+i] S0(3) X R3 X R5 to the ODE (g^) with vector fields Z*,(i?j,s) and Cauchy data 
Rjitj) = R{tj), Yj{tj) = r{tj) and s^{tj) = satisfy: 

1- s^meit,.t,^^] {\\m - Rjimuis) + \\m - r.m^. + - ^t\\ci,iR-^r) < ^72 with §i ■.= 

2. R,{tj+i)l R{t,+li r,(t,+i) = r(<,+i) and {t,+,) = {t,+, - 0, 0, 0, 0)*; 
With these settings, the functions t G [0, T] i-^ i3t G 1^0(1^^)' ^ ^ [0' ^] ^ ^0(3) and f : [0, T] ^ R'^ 
defined by dt ■= ^1, R{t) := Rj{t) and f(<) := rj(t) if t G (j = 1, . . . , fc — 1) are continuous, 

piecewise C^. We extend also every functions on [0,T] by setting (t) := if t G [0,tj[ and 
S''(i) := {tj+i — ij, 0, 0, 0, 0)* if t G]ij+i,r]. They are continuous, piecewise as well. 

It remains to smooth the function -df (and hence also R and f). We can extend the func- 
tions Xl on the whole interval [0,r] by merely setting A^(i) = if t ^ [tj,tj+i]. Then, denoting 
(F^ ) i<,<5 the canonical basis of (R^)*^ and S := Y.'^^i ELi ^l^l ^ i^^f^ we get that (E, f, S) is 

16 



a Caratheodory's solution to the following equation on [0, T]: 



I M = EE^^WT^(^'^)' (5-1) 

\S) j=i«=i 

where T{{R,S) (M.^c^ s^), i?Xf (c^ s^), ) e TflS0(3) x x {H^f. Let denote an- 
alytic approximations of the functions A;^ in L^{[Q,T]) and denote {R,i,S) the corresponding 
analytic solution to system (|5.ip with 5 := (s}, . . . , Sg, s^, . . . , S5, . . . , sj, . . . , Sg) and i9t := i?J + 
Ej=i ELi which is analytic from [0,r] to CqHR^)^- According to Lemma lES - A 

can be made arbitrarily small in AC{[0,T],Cq{'R^)^), providing that the functions A^ are close 
enough to A^ in L^{[0,T]). Notice that does probably not satisfy the self-propelled constraints 
(jl.2p (especially at the times tj, j = 2, . . . , k — 1). It remains to invoke Proposition IB . II and the 
continuity of the input-output function (second point of Proposition II. ip to conclude that there 
exists t e [0,T] n- e £)q(R^) analytic, satisfying (|1.2p and such that 

sup (\\Rit) - mWuis) + ||r(t) - f(t)||R3 + i?t||ci(R3)3) < s/2, 

where (i?, r) : [0,T] 1-^ S0(3) x is the solution to System ((LH) with initial data (i?(0),r(0)) = 
(^(0), r(0)) and control i^. The proof is then complete. 

6. Conclusion. In this paper, we have proved that, for a 3D shape changing body, the ability 
of swimming (i.e. not only moving but tracking any given trajectory) in a vortex free environment 
is generic. The genericity refers to the shape of the body, its density and the basic movements (at 
most five) required for swimming. This result is part of a series of articles [T31 HH HH H] studying 
locomotion in a potential flow and the next step of this study will be to investigate whether this 
controllability result can be extended to a flow with vortices. 

Appendix A. Function spaces. 

Classical function spaces. 

• For any compact set K C R" (n a positive integer), the space C^(K) consists of the 
continuous functions in K endowed with the norm ||w||co(x) sup^^^^ l'"(2;)|- The open 
subset of the positive functions of C^(IC) is denoted C^(K)^ . 

• For any open set f2 C R^ (included = R^), 2?(il) is the space of the smooth (C°°) 
functions, compactly supported in f2. 

• For any open set 17 C R'^ (included 17 = R'^), the set C\(^ is the completion of P(51) for 
the norm ||w||ci(n) •= sup^gQ + II^w(2;)[|r3- When 12 = R^, we get Co(R'^) {u G 
C^(R3) : \u[x)\ and ||Vm(x)||r,3 -J- as ||a;||R3 +00}. 

• The space Cq(R'^)^ is the Banach space of all of the vector fields in R'^ whose every 
component belongs to Cq(R'^). 

• Let E be an open subset or an embedded submanifold of a Banach space and T > 0, 
then AC{[0,T], E) consists in the absolutely continuous functions from [0,T] into E. It is 
endowed with the norm ||w|Uc([o,t],_e) := suptgjg l|ut||£; -l- \\dtUt\\E<it. 

• C™(17, M(fc)) (m an integer) is the Banach space of the functions of class C™ in R^ valued 
in M(fc) (M(fc) stands for the Banach space of the k x k matrices, k a positive integer) and 
compactly supported in 12. 
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• i?™(r2,M(fc)) stands for the connected component containing the zero function of the open 
subset {M e C^{n,M{k)) : det(Id + Af(a;)) t^OVxgR^}. 
Lemma A.l. The set D^{R^) {t? £ C^CR^)^ s.t. Id + is a diffeomorphism of R^} is 
open in Cq(R'^)'^. 

Proof. The mapping i? £ CqHR^)^ := inf ees2 (Id + Vi?(a;), e) ■ e e R (5^ stands for the 

unit 2 dimensional sphere) is weU defined and continuous. For any do G -Dq(R^), we have S^^ > 
and for ah x,y G R'^ and e := {y — x)/\y — x\ the foUowing estimate holds: {y + 'd{y) — x — i?(a;)) ' ^ = 
y — a;| /p (Id + Vi^(a; + te), e) ■ edt > \y — x\S^. We deduce that Id + is one-to-one if d is close 
enough to ^q. Further, still for 'd close enough to i?0; Id + is a local diffeomorphism (according to 
the local inversion Theorem) and hence it is onto. □ 

Definition A. 2. We denote _Dq(R'^) the connected component o/I)q(R^) that contains the 
identically zero function. 

If i9 G Cq(R'^)'^ is such that ||i?||c'i;(R3)3 < li the local inversion Theorem and a fixed point 
argument ensure that Id + is a diffeomorphism so we deduce that Z?q(R^) contains the unit 
bah of Co^(R3)3. 

Sobolev spaces. For any open exterior domain T, the weighted Sobolev space W^{J^) is 
defined by W\T) := {u G V'{T) : u/^/TTW S L^iT), d.^,u e L^{T), i = 1,2,3} (see [3] for 
details). 

Appendix B. Making shape changes allowable. 

Proposition B.l. Let t e [0, T] t-^ e £)q(R^) forj = 1, . . . , +oo, be a sequence of absolutely 
continuous functions (respectively of class C for p = 1, . . . , +cx) or analytic) which converges to 
t e [0,T] € -DolR^) in AC{[0,T], D^R^)). Assume that for some function g G C°(B)+, the 

pair {Qj^^) satisfies (|1.2p . Then, there exists a sequence t € [0,T]i— >t9( G I?o(R"^) (j ^ 1, . . . , +go) 
in AC{[0, T], Dg(R'^)) (respectively of class for p = 1, . . . , -f-oo or analytic) such that (i) for any 
j = 1, . . . , +00, the pair (g, satisfies ([121) and (ii) d'^ i?^ m AC{\0, T],Dl{R^)). 

Proof. Denote m :~ J^gdx and, for every j = l,...,-|-oo, 0^ — Id -I- i?^ and r-'(t) := 
(1/to) Jg gOl dx. Then, for any continuous function t G [0,r] H> R{t) G S0(3), we have /g gR{t)0( dx = 

where Ol := 0^ — r^{t). Let us now determine an absolutely continuous function t G [0,r] t->- 
Wit) G S0(3) such that we have also ^ ^ g[dt{W {t)0{) x m{t)0{]dx = for aU t g]0,T[. In- 
troducing I((9f) := /g£i[|(9fpld ~ 0l ® 0l]dx (an inertia tensor with > 0, so always in- 
vertible) and G{0l,dt0t) := 1{01)~^ j ^ g dt0l x 0(dx, the identity above can be turned into 
R3{t) = W{t)G{0{,dt0l) (0 < t < T). Since S0(3) is compact, this ODE supplemented with 
the Cauchy data i?(0) = Id admits a unique solution defined for all t G [0,T]. The solution has 
the same regularity as t G [0,T] i— > G Dq{R?). Besides, basic estimates allow to prove that 
(i?^r^) ^ (ld,0) in AC([0, T], S0(3) x R3)_and next that i)^ 7?t in ^C([0, T], C(J(R^)^)- No- 
tice that, at this point, we probably have ^ _Dq(R'^). Let now f2 be a large ball containing 
Ute[o,T] 0j (B) and Vl' be an even larger ball containing Q,. Consider x ^ cut-off function such 

that 0<x^liX = lin^^ and x = in R^ \ Q,' . Define 0^ as the flow associated with 
the Cauchy problem X{t,x) = x{x)^t^ + (1 - x{x))dtd\ X(0,x) = 0\=oi^), {x G R^) and 
:= 03 _ Id. Since ^^^^ = we deduce that ^ G Dl{R?) for aU t G [0,r] and the sequence 

t G [0, T] i?^ G Do(R^) complies with the requirements of the Proposition. □ 
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Appendix C. Making vector fields allowable. Let a triplet {g, V) G C"(i?)+ x L'o(R^) x 
(C^(R3)3)" be given such that Jggedx = where = Id + t?. Recall that E = 95. 

Proposition C.l. It is always possible to define new vector fields V* in such a way that (i) 
Vrls = V,|e, (ii) {q, I?, V*) € C(n) with V* := (V^ , . . . , V* ) and (lii) J^gO- V* dx = 0. 

Proof. Arguing like in the proof of Theorem I2.3[ we can easily show that the mapping W G 
Co(-B)'^ H> (Jg gW dx, Jg gOxW dx, Jg gO-W dx) GR'^xR'^xRis onto with infinite dimensional 
kernel. Hence, it is always possible to find a vector field Wi G Cq{B)^ satisfying Jg gWidx = 

— Jg gVidx, Jg gO X Wi dx ~ — Jg gO x Yi dx, Jg gO ■ dx = — Jg gO ■ dx and such 
that {0\b ■ Bk, (Vi + Wi)|b ■ ek, k = 1,2,3} is a free family in C^{Il^f. Wc denote V^; := 
Vi +Wi . We can continue with the same idea: The mapping W G Cq (R"^)"^ i-> (Jg gW dx, Jg gOx 
"W dx, Jg gYl X Wdx) G R'^ x R'^ x R-^ is onto, also with infinite dimensional kernel. Again, 
it is possible to find W2 G Cq{B)'^ satisfying Jg g'W2dx = — Jg gY2 dx, JggO x W2dx = 

— /b C 6' X V2 dx, JggO ■ W2 dx — — Jgg0-Y2 dx and JggYl x W2 dx — — Jg gYl x V2 dx 
and such that {0\b ■ efe, YI\b ■ ek, (V2 + W2)|b • e^, /c = 1,2,3} is free in C^(Il^)^. We can set 
V2 := V2 + W2 and iterate this process to define V3, . . . , V* . □ 

Appendix D. Added mass matrix for particular shaped swimmers. Let B be an open, 
bounded, connected, subset of R'^. Denote U its boundary, n the unit vector to S directed 
toward the interior of i3, := R"^ \ B and consider the 6x6 symmetric matrix of which 
the entries are defined by JjrV^pi ■ 'V'ipj dx (1 < i,j < 6) where the functions ipi (i = 1, . . . ,6) are 
harmonic in J- and satisfy the Neumann boundary conditions dnipi = (e^ x x) ■ n on 17 if i = 1, . . . , 6 
and dnipi = • n if i = 4, 5, 6. 

Proposition D.l. The matrix M.^ is always positive and it is positive definite if and only if 
B has no axis of symmetry. 

Proof. Denote ot := {ai, . . . ,aeY any element in R^ and ip := '^iV'j- Then, we have 

a*Mf a = II V^/'||^3d.T > which proves that M-'' is indeed positive. Let now a be in R^ such 
that oc^M^cx = 0. It means that tp — and hence dnip = oai^i x a;) • n + X]i=4 '^i^i • n = 0. 

Denoting u := '^i^i ^'^^ •= 12^=4 ^i^i • n = 0, this condition reads (u x a; + v) • n = 

on S. The case u = is obviously not possible, otherwise we would have v ■ n(a;) = for all 
X € S s,o let us assume from now on that u 7^ 0. For all a; G S, (u x a: + v) G T^S (the 
tangent space to E at x). Let us set w := (u x v)/|up. Then (u x x + v) = u x (x — w) + Au with 
A (v-u)/|up. In this form, we can explicitly compute the expression of the fiow connecting to the 
ODE X{t,x) = ux (X(<,x)-w) + Au, X{Q,x) = x. We obtain X{t,x) = exp(iu)(x- w) + w + tAu 
for all i G R. Since X{t,x) has to remain on S for alH G R and S is bounded, we deduce that 
Au = 0. For all X G £' and all t G R, the point X{t,x) lies on E. It means that S is globally 
invariant under a rotation whose axis has u as direction vector and goes through w. □ 

Let M. he a. smooth embedded submanifold of a Banach space E and X := (Xi)i<i<„ be set 
of smooth vector fields on AA. 

Lemma D.2. Let ai be in L°°[[0,T]) (i = 1, . . . ,n) and x be a Caratheodory's solution defined 
on the time interval [0,T] to the ODE x{t) = '^^^i cei(t)'X.i(x) (0 < t < T), with Cauchy data 
x(0) — xq M.. Let the functions , (i = 1, . . . , n, fc G Nj be in L^{[0,T]) , such that — > ai 
in L^([0,r]) as k ^ +00. Then, for any sequence of Caratheodory's solutions {xk)keTSi satisfying 
Xk{t) = Q!i'(i)Xi(x) with Cauchy data Xfe(O) = xq, the functions x^ can be continued on the 

whole interval [0,T] for k large enough, x^ — > x uniformly on [0,T] and Xk ^ x in L^{[{),T],M). 

Proof. For any 5 > Q small enough, denote by Kg the compact {x G : ||x — x(i)||£; < 5, t £ 
[0,r]} and denote fc^ > the Lipschitz constant such that X]r=i ll-^i(^) ~ ^i(2/)ll-E < ks\\x — y\\E 
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for all a;, y G ifa and all i = 1, . . . , n. Let 7\/ := max xea-^ ||Xi(a;) jj^; and m := sup fe[o,r] |ai(t)|. 

1 = 1 71 i = l,...,n 

Any function Xk is defined at least on a small interval [0,tfc[ and we can choose tk small enough 
such that Xk{t) £ Ks for all t G [0,tk[- We get the estimate, for all t G [0,tfc[: 

\\x{t) - Xk{t)\\E < \a,{s) - aHs)\ds + mks \\xk{s) - x{s)\\Eds. (D.l) 

Jo Jo 

For any £ > 0, we can choose k large enough such that X]"=i /o Iq^*!*) ~ '^i{^)\ds < £e~™'"'''^/M. 
Applying Gronwall inequality to (jP.ip . we obtain that \\x{t) — Xk{t)\\E < £ for all t G [0,tfc[. We 
deduce first that if £ < 5, the solution Xk can be continued on the whole interval [0, T] and then 
that Xk ^ X uniformly as fc — > +oo. Writing now that: 
pt " pt 

/ \\xk{t) ~ x{t)\\Edt <Y / |a- W -«^Wll|X,(a;fc(s))||£ + |a,(i)|||X,(a;fc(i)) -X,(x(t))|i£dt 
Jo Jo 

" pt pt 

<A/V/ -a,(t)|dt + mfc5 / \\xk{t) ~ x{t)\\Edt, 

»=i -'o Jo 

we get the convergence of the sequence (ifc)fcgN to x in L-^([0, T], A^). □ 
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